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WHOLE-BODY MATHEMATICAL MODEL FOR SIMULATING INTRACRANIAL 

PRESSURE DYNAMICS 

RELATED APPLICATION DATA 

[0001] This application claims the benefit of U.S. Provisional Patent Application 
No. 60/409,551, filed September 10, 2002, entitled "A Mathematical Model for Human 
Intracranial Pressure and Flow Dynamics," and U.S. Provisional Patent Application 
No. 60/416,508, filed October 7, 2002, entitled "A Mathematical Model for Human 
Intracranial Pressure and Flow Dynamics." 

GOVERNMENT LICENSE RIGHTS 

[0002] The U.S. Government may have a paid-up license in this invention and the right in 
limited circumstances to require the patent owner to license others on reasonable terms as 
provided for by the terms of Grant No. DMS-96-26391 awarded by the National Science 
Foundation, and under Grant No. NGT5-401 10 and Cooperative Agreement No. NCC5-581 
awarded by NASA. 

FIELD OF THE INVENTION 

[0003] The present invention relates generally to mathematical models for simulating 
intracranial pressure dynamics and more particularly relates to a whole-body mathematical 
model for simulating intracranial pressure dynamics. 

BACKGROUND 

[0004] Lumped-parameter models represent an attractive method for examining pressure 
dynamics involving complicated human physiology. In a lumped-parameter modeling 
approach, the physiological system is subdivided into a number of linked, interacting subunits 
termed "compartments." In general, each compartment contains a single physical constituent, 
e.g., blood, cerebrospinalfluid (CSF), or tissue and interstitial fluid. However, depending on 
the model's complexity, a given constituent may appear in more than one compartment of the 
model. Dynamics in each compartment is specified by lumped, time-dependent functions 
giving compartmental pressures, while incremental changes in flows and compartmental 
volumes are obtained by associating resistance and compliance parameters with adjacent 
compartments. In particular, interaction between adjacent subunits is assumed to take place 
at the interfaces of the model's compartments. 



[0005] With few exceptions, previous models of this type have adopted restrictions known 
as the "Kellie-Monro Doctrine" to reduce complexity. The Kellie-Monro framework 
considers the intracranial system to be completely enclosed within the intracranial vault, 
which is assumed to be rigid. A specified inflow of blood to the intracranial arteries provides 
a forcing for the system, and outflow from the jugular bulb is assumed to instantaneously 
equate to this inflow. These restrictions yield a closed system with constant total volume. 
Strictly intracranial models have produced a number of important results that illuminate the 
mechanisms of intracranial pressure adjustments in situations involving both normal and 
pathophysiology. However, the ability of these closed-system models to incorporate the 
influence of important extracranial factors on intracranial pressure dynamics is clearly 
limited. For example, the important buffering effects of the spinal CSF space on intracranial 
pressure cannot be directly included. From a mathematical point of view, the constant 
volume constraint also produces an over determined system of equations that requires special 
handling to avoid singular behavior. 

[0006] Lumped-parameter compartmental models of the present type have a long history, 
dating to the earliest such model of the intracranial system formulated by Monro in 1783. 
This first model was bi-compartmental, considering incompressible brain matter and blood as 
its two constituents. In the work of Monro's student Kellie 40 years later, the vascular 
compartment was further subdivided into arterial and venous blood to produce a three- 
compartment model. Since the pioneering work of Monroe and Kellie, increasingly more 
complex models of the intracranial compartment have been posited in order to more 
realistically describe the relationship between intracranial pressures and volumes. There has 
been a steady increase in the number of fluid compartments, the introduction of a separate 
cerebrospinal fluid compartment, the inclusion of cardiovascular input, and a relaxation of the 
treatment of system constituents as incompressible fluids and matter. As noted in prior art 
studies, the intracranial system involves a number of subsystems, which interact through 
complex mechanisms, and the classical piecewise approach, which often considers CSF 
dynamics separately, is not suited to studying multiple parameter changes and the effects of 
interconnected subsystems on each other. By contrast, lumped-parameter models of the 
intracranial system are capable of including and linking different subsystems, so that such 
interactions can be examined. 
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[0007] When considering lumped parameter models, it is important to realize that a 
compartment does not necessarily correspond to a precise physical location in the body. For 
example, with a single CSF compartment in a model, CSF in the ventricles cannot be 
distinguished from CSF in the subarachnoid and spinal spaces. This is one of the main 
limitations of the lumped-parameter approach. Additional spatial resolution can be realized 
only by subdividing the physical system into a larger number of compartments based on 
spatial considerations. For example, distinct ventricular CSF and extra- ventricular CSF 
compartments may be included as opposed to a single lumped CSF compartment. In 
principle, the entire body could be finely subdivided in this manner into separate 
compartments to provide the desired degree of spatial resolution. However, clearly this 
subdivision process cannot be carried to an extreme as the resulting system of linked 
governing equations will rapidly become too large for practical analysis and solution. 

[0008] Despite their evolving complexity, two common features characterize most earlier 
lumped parameter models for pressure dynamics in the intracranial system. The first 
common feature is an assumption that all resistance and compliance parameters can be 
represented by constants. This leads to a linear system of governing differential equations. 
The second common feature is adoption of the "Kellie-Monro Doctrine," which assumes that 
the intracranial system can be confined within the cranial vault. By requiring that inflow to 
the intracranial arteries equals outflow from the jugular bulb, this assumption produces a 
closed system that conserves total intracranial volume. 

[0009] As mentioned above, when the intracranial space is treated as a closed volume- 
conserving system contained within the (nearly) rigid cranial vault, important mechanisms for 
the influence of extracranial physiology on intracranial pressure dynamics cannot be included 
in the resulting models. For example, the ability of the spinal portion of CSF space to buffer 
fluctuations of intracranial CSF pressures cannot be directly introduced under the Kellie- 
Monro Doctrine. At least two prior art mathematical models of intracranial pressure 
dynamics include aspects of extracranial physiology. The first model includes a compliance 
for CSF storage within the lumbar channel. The first model contains three compliances, four 
resistances, and involves differential equations based on a hydrodynamic model for the 
physical system and its electrical circuit equivalent. The first model allows the dynamic 
relationship between cerebral perfusion pressure, intracranial pressure, and cerebral blood 
flow in various states of autoregulation to be studied. Use of the first model in conjunction 
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with clinical data has determined which indices can be derived using transcranial Doppler 
ultrasonography and which trends of intracranial pressure and blood pressure are useful in 
clinical tests of autoregulatory reserve. However, despite not strictly abiding by the Kellie- 
Monro Doctrine, the first model falls far short of being a whole-body model for intracranial 
pressure dynamics. While the first model can include direct buffering effects of CSF within 
the lumbar channel, it does not contain other important extracranial physiology. For 
example, the first model fails to include representations for the extracranial physiology of the 
cardiovascular system and the tissues and extracellular fluids of both the central and lower 
body, including cardiovascular autoregulation, colloid osmotic pressure, and a lymphatic 
system. 

[0010] The second model uses a single ground compartment to represent the portion of the 
body below the clavicles. The second model contains three resistances, including a resistance 
between the intracranial veins and the rest-of-body compartments. Thus, outflow from the 
intracranial system depends on pressure differences and need not instantaneously equal the 
specified inflow. The second model also contains three compliances, including a compliance 
between the CSF and rest-of-body compartment that represents the ability of the spinal 
portion of CSF space to expand and buffer CSF pressures by compressing the large veins that 
surround the spinal theca. Two of the three model compliances are pressure dependent. The 
associated functional forms for the pressure-dependent compliances determined in the second 
model are used in the present work. The second model still falls short of being a whole-body 
model for intracranial pressure dynamics. The second model merely lumps extracranial 
physiology into a single compartment, rather than finely subdividing it into a plurality of 
compartments. In addition, the second model fails to include the external environment as an 
implicit compartment. 

SUMMARY OF THE INVENTION 

[0011] One aspect of the present invention is a method of modeling pressure dynamics of a 
body's intracranial system including the following steps: dividing the body into a plurality of 
compartments and a heart pump, each of the plurality of compartments representing a portion 
of the body, the heart pump interacting with at least one of the plurality of compartments; 
deriving a plurality of differential flow equations, each of the plurality of differential flow 
equations corresponding to one of the plurality of compartments; and solving the plurality of 
differential flow equations. 
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[0012] Another aspect of the present invention is a system for modeling an intracranial 
system including a body compartment module adapted to divide the body into a plurality of 
compartments and a heart pump, each of the plurality of compartments representing a portion 
of the body, a flow equation module adapted to derive a plurality of differential flow 
equations, each of the plurality of differential flow equations corresponding to one of the 
plurality of compartments, and an equation solver module adapted to solve the plurality of 
differential flow equations. 

[0013] Still another aspect of the present invention is a method of modeling pressure 
dynamics of an intracranial system that further accounts for cerebrovascular autoregulation 
by a sympathetic nervous system. 

[0014] Yet another aspect of the present invention is a method of modeling pressure 
dynamics of a body's intracranial system including the following steps: dividing the body 
into a plurality of compartments and a heart pump, each of the plurality of compartments 
representing a portion of the body, the heart pump interacting with at least one of the plurality 
of compartments, wherein a plurality of the plurality of compartments are vascular and a 
plurality of the plurality of compartments are non-vascular, the vascular compartments 
including at least one of the intracranial arteries, intracranial capillaries, choroids plexus 
capillaries, venous sinus jugular veins, intracranial veins, central arteries, central capillaries, 
central veins, extra-ventricular CSF, lower arteries, lower capillaries, and lower veins, the 
non-vascular compartments including at least one of lower tissue, brain, ventricular CSF, and 
extra- ventricular CSF; deriving a plurality of differential flow equations, each of the plurality 
of differential flow equations corresponding to one of the plurality of compartments; and 
solving the plurality of differential flow equations. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] For the purpose of illustrating the invention, the drawings show a form of the 
invention that is presently preferred. However, it should be understood that the present 
invention is not limited to the precise arrangements and instrumentalities shown in the 
drawings, wherein: 

FIG. 1 is a simplified schematic of one embodiment of a multi-compartmental model 
for determining intracranial pressure dynamics; 
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FIG. 2 is a block diagram of the various compartments of a model according to 

FIG. 1; 

FIG. 3 is the intracranial arteries compartment of the model according to FIG. 1; 
FIG. 4 is a chart of interstitial fluid pressure vs. relative lymph flow; 
FIG. 5 is a chart of percent reduction in Q AC vs. P A ; 

FIG. 6 is a simplified schematic of a computer system including one embodiment of 
the present invention; 

FIG. 7A is a chart of pulsatile cardiac forcing (Qhi) over time; 

FIG. 7B is a chart of the 120/80 central artery pressure response over time (solid) and 
the i 00/65 intracranial artery pressure response (dashed); 

FIG. 8A is a chart of the predicted CSF pressure vs. volume; 

FIG. 8B is a chart of conductance to CSF output according to one embodiment of the 
present invention; 

FIG. 9 is a chart of the percentage of Q AC drop due to a positional change with all 
SNS reflexes activated (solid) and without any SNS reflexes included (dashed); 

FIG. 10 is a chart of percentage decrease of P \ (solid) and of Q H i (dashed) as a 
function of percentage blood loss due to hemorrhage; and 

FIG. 1 1 is a chart of the predicted percentage of Q A c drop vs. percentage of blood loss 
due to hemorrhage with both SNSz reflexes activated (solid) and no SNSz reflexes included 
(dashed). 
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DETAILED DESCRIPTION OF THE INVENTION 

[0016] The present invention is directed to a whole-body mathematical model for 
simulating intracranial pressure dynamics. Examples of the disclosed invention are depicted 
in FIGS. 1-11, although it should be understood that the present invention is not limited to 
this (or any other) particular embodiment, but rather is intended to cover models and systems 
that fairly fall within the broad scope of the appended claims. 

[0017] As described in U.S. Provisional Patent Serial No. 60/409,55 1 , which is 
incorporated by reference as if disclosed herein in its entirety, the model of the present 
invention revokes the Kellie-Monro Doctrine by consistently embedding the intracranial 
system within whole-body physiology. In one embodiment, the model of the present 
invention includes 17 interacting compartments, of which nine lie entirely outside of the 
intracranial vault. Two distinct compartments are defined to distinguish ventricular from 
extraventricular CSF. The vasculature of the intracranial system within the cranial vault is 
also subdivided into five compartments representing fluid in the intracranial arteries, 
capillaries, choroid plexus, veins, and venous sinus. The body's extracranial systemic 
vasculature is divided into six compartments representing the arteries, capillaries, and veins 
of the central body and the lower body. Tissue and the associated interstitial fluid in the 
intracranial and lower regions are divided into two compartments. A composite compartment 
involving the tissues, organs, and pulmonary circulation of the central body and an implicit 
compartment representing the external environment complete the model. Since the time- 
dependent compartmental pressure functions are obtained from physical pressures through a 
"lumping" procedure that involves space-averaging over the physical extent of a 
compartment, the subdivision of physical constituents into distinct spatially limited 
compartments is necessary to provide spatial resolution in this modeling approach. 

[0018] In addition to allowing direct flows (e.g. arteries to capillaries to veins) between 
adjacent compartments, the present work includes the extracranial transfer of fluid between 
capillaries and tissue by filtration. An extracranial lymphatic system is also included in the 
model. Components of the model allow regulation of systemic vascular pressures by the 
sympathetic nervous system, and, at less than extreme (high or low) pressures, autoregulation 
mechanisms provide constant blood flow to the cerebrovascular capillaries and the choroid 
plexus as well as constant production of CSF in the choroid plexus. Fluid intake, renal output 
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of fluid, and adjustment of body volume in response to changes in ambient environmental 
pressure are allowed. A realistic representation for cardiac uptake and cardiac output 
provides the forcing for this system. 

[0019] The form of the present model is a system of governing differential equations for the 
fully time-dependent compartmental pressure functions. A component of the present work is 
determining appropriate forms for the non-constant resistance and compliance parameters in 
the model, which may be functions of both pressures and time. Calibration of physically 
realistic scale values for parameters and flows is also a step in the simulation process. 

[0020] Turning to the drawings, wherein like reference numerals refer to like elements, 
FIG. 1 illustrates one embodiment of a mathematical model 10 according to the present 
invention. In mathematical model 10, the human body is divided into 16 distinct 
compartments, lettered A, B, C, D, F, G, I, J, O, P, S, T, V, X, Y, and Z. Compartment A 
represents the intracranial arteries. Compartment B represents the brain. Compartment C 
represents the intracranial capillaries. Compartment D represents the lower capillaries. 
Compartment F represents the ventricular CSF. Compartment G represents the lower tissue. 
Compartment I represents the central arteries. Compartment J represents the central 
capillaries. Compartment O represents the central veins. Compartment P represents the 
choroid plexus capillaries. Compartment S represents the venous sinus jugular veins. 
Compartment T represents the extra-ventricular CSF. Compartment V represents the 
intracranial veins. Compartment X represents the lower veins. Compartment Y represents 
the rest of the body. Compartment Z represents the lower arteries. Compartment M 
represents the external atmosphere. Compartment H L and H R represent the left and right 
heart chambers. The physical constituents in subunits of the present model are blood, CSF, 
and tissue and interstitial fluid. With a single exception, i.e., compartment Y, each 
compartment is composed of a single constituent. Spatial resolution is obtained by first 
dividing the body into intracranial and extracranial components. To help delineate the 
intracranial system, a thick line 12 in FIG. 1 represents a cranial vault 14, which differentiates 
intracranial from extra-cranial compartments. Compartments I, Z, Y, J, D, G, O, and X lie 
completely outside of cranial vault 14, and compartment T includes both intracranial and 
extracranial physiology. 
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[0021] As illustrated in FIG. 2, the majority of the compartments in the current model are 
vascular. The 1 1 vascular compartments can be subdivided spatially into three groups: 
intracranial compartments A, C 9 P, V, and S; central compartments I, J, and O, and lower 
compartments Z, D, and X. The "lower" compartments represent the region below the pelvis. 
The "central" compartments in the model represent the region between the lower body and 
the clavicles and also include extracranial body components above the clavicles but outside 
of the cranial wall. The vascular system in each of these three regions, i.e., intracranial, 
central, and lower, is divided into separate artery, capillary, and vein compartments. 
However, in the intracranial space, the choroid plexus capillaries are placed in a separate 
compartment P from the rest of the intracranial capillary bed. This allows these capillaries to 
autoregulate to maintain the production of CSF when ventricular CSF pressures rise or fall. 
The venus-sinus veins where CSF absorption occurs through the arachnoid villa and the 
jugular veins are also placed in a compartment separate from the remainder of the intracranial 
veins, i.e., compartment S. There are four strictly non-vascular model subunits. Two of 
these compartments represent tissue matter and interstitial fluid, i.e., lower compartment G 
and brain compartment B, and the other two represent CSF compartments, i.e., ventricular F 
and extra-ventricular T. Compartment T is both an intracranial and a central compartment. 
Compartment T, which contains extra-ventricular CSF and lies partially outside of the cranial 
vault, includes both the subarachnoid and spinal portions of CSF space. Compartment T 
serves as a bridging compartment and explicitly allows for buffering of CSF pressures by the 
spinal theca. Two of the three regions contain a separate compartment for the tissue and 
interstitial fluid. The exception is the central region where the tissue matter, organs (except 
for the heart), interstitial fluid and pulmonary circulation are lumped into a composite rest-of- 
body compartment Y. The central region contains an explicit representation for the heart 
pump having a left chamber pump H L , and a right chamber pump H R , and a realistic pulsatile 
cardiac output provides the major forcing in the current model. External environment 
compartment M is also an implicit 17-th subunit in the model. No attempt has been made to 
depict relative volumes in FIG. 1, and hence the relative sizes of these compartments in the 
figure do not reflect relative volumes. 

[0022] The pressure dynamics of the intracranial system are governed by a system of 
differential equations within mathematical model 10. Four basic assumptions lead to these 
equations: 
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(1) all fluids are considered incompressible and isothermal; 



(2) pressure driven flows are laminar and related to pressure differences by 

P — P 

Qu = J ^ L = z u i p i - p j)= z u p u> ( i ) 

where Q 0 is the flow from compartment i into compartment j, P, and Pj are the spatially- 
averaged pressures of compartments i and j respectively, R 0 is the lumped resistance, Z 0 is the 
fluidity (inverse of Ry) 9 the pressure difference P 0 = P, -Pj , and R 0 = 

(3) in the case of fluid filtration from the capillaries into the interstitial space, the 
flow is governed by the Starling-Landis Equation, i.e., 

Filtration = K ct ((P c - P t ) - a ct (tt c - n t )) = K ct (P ct - <j ct n ct ) , (2) 

where P c is the capillary pressure, P t is the interstitial fluid pressure, n c is the blood plasma 
colloid osmotic pressure, % t is the interstitial fluid colloid osmotic pressure, K ct is the filtration 
coefficient, and o cl is the capillary membrane reflection coefficient. The notation for pressure 
difference has been extended to osmotic pressure differences by defining tt c , = n c - n t \ and 

(4) the deformation of the membrane between adjacent compartments is a 
function of the change in pressure difference between these compartments, i.e., 

dt ^ dt ~ Lij ^r (3) 

where V 0 denotes the instantaneous volume of the 'cup' formed in the membrane at the 
interface of compartments / and j, C 0 denotes the compliance between these two 
compartments, and C 0 = C Jt . 

[0023] As follows, the system of governing equations is obtained by invoking the 
conservation law 
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flow rate in - flow rate out = rate of volume change 



(4) 



in each compartment. Referring now to FIG. 3 5 compartment I from FIG. 1 is illustrated. 
Compartment I represents the central arteries. The flow rate into compartment I is 
represented by Q H i from left heart chamber pump H L . The flow rate out of compartment I is 
the sum of the flows associated with lumped resistances Ri Aj Ru, and Riz, as represented by 
equation (1). The rate of volume change of compartment I is equal to the sum of the terms 
associated with the compliances Crr and Ciy, as represented by equation (3). 

[0024] Specifically, the pressure dynamics of compartment I are governed by the following 
equation: 



[0025] Note that in equation (5), the subscript "IY" has been reversed thereby causing the 
associated term to be negative. Also, in both FIGS. 1 and 3, a filled arrow indicates a one- 
way flow and a hollow arrow indicates a pressure dependent resistance. In terms of pressure 
differences, the relation of equation (4) yields the following additional 12 differential 
equations and three scalar equations for the remaining 15 body compartments: 

- Central Capillary Compartment (J): 



Qw -( z ij p u + z iz p iz + z ja p /a) z:zC 



IT 



dP JT dVyj 
dt dt 



(5) 



z u p u ~{ k jy( p jy - <y jyx jy)+ z jo p jo) = 0 s 



(6) 



- Rest of Body Compartment (Y): 



Qw + k jy( p jy jyXjy)-(Qyo + Qk) = c ym 



d p YM , dV Y1 s dV YO m 



;(7) 



dt dt dt 



- Central Venous Compartment (O): 



z jo p jo + z xo p xo + z so p so + z to p to + Qgo + Qyo - Qoh = C TO 



dPgj dVyp 
dt dt 



; (8) 
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- Lower Artery Compartment (Z): 



Z IZ P IZ ~ Z ZD P ZD ^ C ZG—~ I 

at 



- Lower Capillary Compartment (D): 



Z ZD P ZD ~ ( Z DX P DX + K DG ( P DG " ° ' DG n ' DG )) = 0 '•> 



- Lower Tissue Compartment (G): 



&DG ( P DG - ct dg tt dg )-Q GO = C zg ^f- + C GX 



- Lower Venous Compartment (X): 



Z DX P DX ~ Z X0 P XO = C GX I 

at 



- Intracranial Artery Compartment (A): 



Z, a P,a - {ZacPac + ZapPap ) = C AB + C AF ^f- 

at dt 



- Intracranial Capillary Compartment (C): 



Z AC P AC ~( Z CB P CB + Z CV P Cv) = ®\ 



- Choroid Plexus Compartment (P): 



Z AP P AP ~{ Z PF P PF + Z PV P Pv) z=C PF ~~ \ 

at 
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- Intracranial Veins Compartment (V): 



dP HP 

Z C yP C y +Zp V Ppy +Z B yP B y - Zy S Py S = C B y ^ + C FV ^fE. J (16) 

at at 



- Venous Sinus - Jugular Veins Compartment (S): 



Z VS P VS + ^75^ ~ Z SO P SO =CTS ' (17) 

at 



- Ventricular CSF Compartment (F): 

dP F j „ dP FP „ dP™ _ dP 



"* «f a/ dr 



(18) 



- Extra-Ventricular CSF Compartment (T): 

(19) 



ZftPft + Z 57 -/> flr - (Z rs P ra + Z ro /> ro ) = C J3 ^f- + C BT ^ + C T0 ^L + C IT ^P-- and 



- Brain Compartment (B): 



„ dP RV „ rfP DC . dP Br 



Wca + Z ra P ra - (Z 5K /> BK + Z BT P BT ) = C AB ^M. + C BV ^f- + C FB ^ + C BT . 

at dt . dt dt 

(20) 

[0026] The terms dVy/dt and dV Y c/dt in the conservation equations for compartments Y, I, 
and O have been left intact in equations (5), (7), and (8). These volume changes include 
components that reflect the regulation of arterial pressure by the sympathetic nervous system. 
They are considered further below. The terms involving dP^/dt in the conservation equations 
for 6 compartments Y, and G (equations (7) and (11)) denote a volume change for which 
there is no compensation through an equal but opposite volume change in an adjacent body 
compartment. Specifically, these terms reflect volume increases (or decreases) into the 
ambient environment, which is considered infinitely large and unaffected by pressure changes 
within the body. Ambient environment compartment M is similar to the ground in an 
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electrical circuit model. It will be maintained at a constant pressure value in exactly the same 
way that the ground voltage in a circuit model is given a prescribed value relative to which all 
other values are gauged. Accordingly, the ambient environmental pressure Pm is fixed here at 
a constant reference value of zero mmHg. If this is not appropriate for a situation being 
studied, or if the ambient pressure changes with time, then P M must be considered as a 
forcing term in the system of equations. 

[0027] The set of governing equations may be summed to derive the following constraint 
regarding compliances between the body and the external environment: 



(Qw -Qk)+ (Qhj ~Qoh)= C Y m + C GM *<jj*L . (21) 

at at 



[0028] If cardiac output equals cardiac uptake (Q HI = Qoh) and the atmospheric pressure is 
constant, equation (21) simplifies to 



dPy 

Qw -Qk = C YM ~~T~ + C GM ~~j • (22) 

at at 



[0029] If, in addition, Q w = Q K = 0 or fluid intake equals fluid discharge (Q w = QkX then 
equation (22) implies as expected that the net volume change of the entire body system must 
be zero. Furthermore, if Qw > Q K > as will initially be the case when a glass of water is 
consumed, the net volume of the body will increase. By equation (22), one of the internal 
compartmental pressures: Pq, or most likely Py, must now increase. However, if as will be 
the case, Cym is large, the resultant pressure increase in compartment Y due to the volume 
input from the ambient environment will be small. 

[0030] The fluidities in equation (1) must be derived, the filtration and reflection 
coefficients in equation (2) must be calibrated, and the compliances in equation (3) must be 
calibrated before model 10 5 s governing differential equations may be used in simulations. 

[0031] Each constant fluidity Zy between arbitrary adjacent compartments / and j may be 
calculated by substituting the mean pressures and the mean flow rates into equation (1) and 
solving for Zy, giving 
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[0032] Therefore, once the mean flows and pressures are estimated from physical data, the 
associated constant fluidities may be calculated from equation (23). Some model fluidities 
are pressure dependent. A discussion of appropriate expressions for pressure dependent is 
provided below where CSF, cerebrovascular, sympathetic nervous system, and cardiac 
autoregulation mechanisms are modeled. However, even for a non-constant fluidity, a mean 
scale value may still be calculated from equation (23). In the calibrations that follow, mean 
flows and pressures reflect physiological values for an average human in a lying down 
(supine) position. 

[0033] To estimate mean flows, it will be assumed that compartmental volumes remain 
constant in the mean state. Compartmental pressures are certainly pulsatile, so the system 
will have a mean state, but not a steady state. However, since volume changes are related to 
pressure differences between adjacent compartments, if pressures in adjacent compartments 
move in tandem in the mean state, volumes will remain unchanged. Therefore, for a given 
compartment, if all but one of the mean flows are known, the final mean flow may be 
determined from maintaining constant volume. Further, once mean flows into and out of a 
given compartment are determined, these values provide data for flows in the adjacent 
compartments. While many flows can be estimated from data available in published 
literature, most mean flow calibrations must make use of the constant volume assumption. 

[0034] As a starting point for the calibration of mean flows, consider the percentages of 
cardiac output Qhi that exit the central arteries into the three peripheral artery systems: A, I, 
or Z. These percentages, in decimal form, are given below. 

Pia = . 15, = % of cardiac output into intracranial region (24) 

p iz = .25, .35, = % of cardiac output into lower region (25) 
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Pij = 1 - {pia + piz), - % of cardiac output into rest of body (26) 



[0035] Additional helpful percentages and ratios include 



p P f - .70, = % of CSF formation from the choroid plexus (27) 



?, = 250 5 = Qac'Qap (28) 



p ts = .80, = % of CSF drained into venous sinus (29) 



[0036] Literature values give the following mean flows in ml/min: 



Q H} = 5000, 6600, 6900, Cardiac output (30) 



Q CF = 2, Total capillary filtration (with 2/3 from liver and intestines) (31) 



QjY ~ Qcf 



%1 El " 



^3 ?>Pij+Pi- 



, Total capillary filtration into central tissue (32) 



) 



1 & 

Qdg ~ " Qcf > Total capillary filtration into lower tissue (33) 

3 Pij + Piz 



Q F = .35, Total CSF formation (34) 



Qbv ~ -001 > Imbalance of diffusion on venule side of cerebrovasculature (35) 



Qfb = - 044 ' Flow of CSF through the Virkow-Robins Space (36) 

[0037] All of the remaining mean flows in the model may now be calibrated based on 
equations (24) through (36) by invoking the constant volume assumption in each 
compartment during the mean state. The calibration procedure yields the mean flows: 
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Central Body Mean Flows 



Qoh - Qhi fr° m cardiac output = cardiac input (37) 

QlA = PiaQHI from equation (24) (38) 

QlZ = PhQhi from equation (25) (39) 

Qu = Qhi ~ {Qia + Qiz ) from constant volume in I (40) 

Qjo ~ Qu ~Qjy f rom constant volume in J (41) 

Qyo = Qjy f rom constant volume in Y (42) 
Lower Body Mean Flows 

Qzd ~QlZ fr° m constant volume in Z (43) 

Qgo -Qdg f rom constant volume in G (44) 

Qdx = Qzd ~ Qdg fr° m constant volume in D (45) 

Qxo ~ Qdx f rom constant volume in X (46) 
Intracranial Mean Flows 

Qap = Qia /(k + 1) from constant volume in A and equation (28) (47) 

Qac = AQap from equation (28) (48) 
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Qpf = PpfQf f rorn equation (27) (49) 



2cs = (l " P p f )Qf + 0£K from g C5 - fi^ = (l - p p/ )a 



(50) 



Qcv = £?/ic ~ Qcb fr° m constant volume in C (51) 



Qpv - Qap - Qpf fr° m constant volume in P (52) 



Qvs = Qcv + Qpv + Qbv fr° m constant volume in V (53) 



Qbt ~ Qcb + Qfb ~ Qbv fr° m constant volume in B (54) 



Qft ~ Qpf ~ Qfb fr° m constant volume in F (55) 



^TS = PisQf from equation (29) (56) 
Qto = 0 ~ )2f fr° m constant volume in T (57) 
= Qvs + 2rs fr° m constant volume in S (58) 

[0038] Equation (49) states that 70 % of the CSF formation comes from the choroid plexus 
and equation (50) states that the remaining CSF formation comes as a filtrate from the 10 
capillaries through the brain via the Virchow-Robins system minus the amount reabsorbed at 
the venule side of the capillaries {Qbv). Notice that, as should be the case, 

Qbt + Qft - Qf • 

[0039] In addition to allowing direct flows (e.g. arteries to capillaries to veins) between 
adjacent compartments, the present model includes the transfer of fluid between capillaries 
and tissue by filtration. These flows are governed by the Starling-Landis equation and are 
driven by differences between the colloid osmotic pressures of the blood plasma in the 
capillaries and the interstitial fluid as well as by the usual compartmental pressure differences 
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between capillaries and tissue. Filtration mechanisms are included between the capillaries 
and tissue compartments in the central and lower portions of the body in the present model. 
In the intracranial region, significant colloid osmotic pressure differences do occur between 
the intracranial capillary and tissues. However, the endothelial cells that make up the 
intracranial capillary wall are so tightly joined that not even water molecules can usually pass 
between them. Thus, under normal conditions, colloid osmotic pressure forces in the 
intracranial region are irrelevant. In the case of highly elevated capillary pressure it is 
possible for the intracranial capillary wall to expand enough to allow water molecules to pass 
between the endothelial cells, and at this point both colloid osmotic pressures and regular 
pressures start governing filtration rates. At these elevated pressures, the description of 
intracranial capillary filtration as a function of capillary pressure will clearly be nonlinear. 
However, to simplify the equation for intracranial filtration, this relationship may still be 
linearly approximated with the slope defined by Z c & Z B v, and Z FB . 

[0040] There are two locations in the model where fluid flow is governed by the Starling- 
Landis equation (2). These flows in the central and lower body are 

Qjy = Kjy {{Pj -Py)- <?jy (*j ~ *Y )) and (59) 

Qdg = K DG ((P D - P G )-a DG (tt d -tt g )), (60) 

denoting the flow from the capillaries into the tissue subunits of the central and lower body 
regions, respectively. Solving these relations for the filtration coefficients in the mean state 
implies 



KjY ={{Pj-Py) Q Zv{*j-*v)) (61) 



K DG = W = v Qd ° ^ (62) 

{{ p D -Pa)- ^dg [*d - n c )) 
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[0041] The mean flows on the right hand side of these equations have been calculated 
above. Thus, to determine the filtration coefficients K JY and K D g, it is only necessary to 
calibrate scale values for the mean colloid osmotic pressures and the reflection coefficients. 



[0042] Mean interstitial fluid colloid osmotic pressure is given by 



7i y = 7iQ = 8 mmHg 



(63) 



and blood plasma colloid osmotic pressure by 



7i j - tt d = 28 mmHg. 



(64) 



[0043] Notice that these values are invariant with respect to central or lower regions. This 
is not the case, however, with the reflection coefficients a JY and o D g> The reflection 
coefficient 1 1 in the legs is estimated to be approximately 0.9 while the coefficient of the 
upper body is less than this value. This is reflected by the assignments 



The filtration coefficients K JY and K DG in equations (61) and (62) are now readily determined. 

[0044] Calibration of scale values for all resistance and compliance parameters from 
available physical data and other relationships must be accomplished before the model's 
governing differential equations can be used in simulations. A key step in calibrating model 
compliances is determining the distensibilities of the extracranial compartments of the model. 
In particular, each extracranial compartment has an associated total compliance equal to the 
product of its associated volume and distensibility. It should be noted that calculating 
compliances by this technique yields a total central artery compliance Q= 1.529ml/mmHg, 
which is within 5 percent of the data value of 1 .445 ml/mmHg for total arterial compliance 
measured in the ascending aorta by Chemla et al. The central venous compliance calculated 
by these methods yields C v = 50 ml/mmHg while known systemic venous pressure volume 



(65) 



odg = 9 



(66) 
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curves suggests a total venous compliance of 50 ml/mmHg. Finally, Noordergraaf states that 
the systemic arterial compliance is between 1 and 2 ml/mmHg and systemic venous is 
between 50 and 200 ml/mmHg. 

[0045] In this section, the various compliances in equation (3) that relate volume 
adjustments to pressure differences will be calculated. In the intracranial region, 
compartmental volume increases are restricted by the rigid cranial wall. Consequently, 
compartmental compliances must be pressure dependent and diminish to zero as pressure 
differences become large. On the other hand, in extracranial regions, to lowest order 
compliances may consistently be considered constant, approximating a linear relationship 
between pressure differences and volume adjustments. The present intracranial pressure- 
dependent compliances are extensions of those derived in the prior art, while the constant 
extra-cranial compliances will be derived from estimations of the volume and distensibility of 
each compartment. 

[0046] In a simplified four-compartment whole-body model for CSF dynamics described in 
the prior art, there are only two pressure difference dependent compliances. They allow 
volume adjustments between the CSF and arterial blood and between the CSF and venous 
blood. These compliances have the general form described by the relation 




(67) 



where Py = Pi - Pj , and the subscripts / and j take the values CL.(for arterial/CSF 

compliance) and c\ (for CSF/venous compliance). For both pairs of index values, 
Coefficients and parameters in equation (67) are given by 




6.5333, 



/> = 0.633431 



atfi = 0.604229 



(68) 




1.82745, 



r af = 0.817102 



Oaf- 0.869393 



(69) 
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[0047] This four-compartment model contains an additional constant compliance between 
the CSF and the rest of the body (g), which is approximated as 



C* g = 0.13333. (70) 



[0048] This feature represents the interface of extra-cranial CSF in the spinal theca with the 
rest of the body. It also acts as a background compliance so that the total CSF compliance 
can never reach machine zero in simulations. 

[0049] In the present 16 compartment model, the division of the body's cerebrospinalfluid 
(CSF) space is considerably more refined. Hence, the three CSF compliances in the simple 
four compartment model of the prior art must be appropriately apportioned among the 
compliances of the present model. Three decimal percentages may be introduced to describe 

this allocation of : 



pjv = .164 percentage of allocated to C FV (71) 



p to = .214 percentage of allocated to C T o (72) 



p ts = .622 percentage of allocated to Crs (73) 



[0050] Equations (71 - 73) reflect the fact that total CSF volume is approximately 140 ml, 
23 of which is found in the ventricles, 30 in the spinal cord subarachnoid space (theca) and 
the remainder in the cerebral cisterns and subarachnoid space. Thus, if the distensibility of 
the membrane is similar in these three different components of CSF, then Cp is 23/140 = .164 

of C^as indicated by p/, in equation (71). This same technique yields the values in 

equations (72) and (73). The compliance C PF is excluded in this division of as the 

choroid plexus capillaries are known to dilate and constrict in order to maintain a constant 
pressure difference between these capillaries and ventricular CSF. This maintains a constant 
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generation rate of CSF from the choroid plexus at normal pressure values. Consequently, a 
value for the compliance C PF is irrelevant in the absence of a pressure difference change. 



[0051] The ratios that lead to equations (71) - (73) also imply the following percentages 
describing the allocation of : 



p a /= .786 percentage of allocated to C AF 



(74) 



pu = .214 percentage of allocated to C !T 



(75) 



[0052] The background compliance C* g is now divided between venous and arterial 
interfaces based on systemic venous volume being four times that of arterial volume, giving 



(-'venous • 8 C 



4 

fg 



(76) 



C arterial «2 C 



fg 



(77) 



[0053] The CSF-related compliances may now be calculated based on the above 
percentages as: 



C FV (Pfv) = .95 p fi (C^ (P FV ) + C venous ) Cfv = 0.557868 ml/mmHg 



C FB (Pfb) 



05 Pfi ( C^, (Pfb) + C venom ) C fb = 0.036255 ml/mmHg 



Cts(Pts) 



0.95 p ls ( (P ra ) + C venous ) Cts = 1 -27626 ml/mmHg 



Cbt(Pbt) = 0.05 p, 5 (C*, (P Br ) + C venous ) C B T= 0.137057 ml/mmHg 
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C T o (Pro) =Pto (Cj, {P TO ) + C venous ) Cm = 0.200936 ml/mmHg 

Caf(Paf) = Paf(C* f (P AF ) + Caneriai) C af = 0.026 1 999 ml/mmHg 

C/r f?VJ =p it (C* f (P /T ) + Carteriai) C it = 0.00571427 ml/mmHg 

where again, the pressure difference Pi - Pj is denoted P 0 and the scale value Cy = 

~^ /) = C^(Pg r). The compliance C B r, introduced to incorporate volume adjustments 

between the brain and subarachnoid CSF, is allocated 5% of C^. This states that of the bulk 
intracranial compliance in compartment T, 95% is allocated to the interface with the highly 
distensible venous sinus veins and only 5% to the interface with the less distensible brain 
tissue. A similar allocation is made for the bulk intracranial compliance of the ventricular 
CSF compartment F with respect to the intracranial veins and the brain tissue. 

[0054] Finally, brain/vascular compliances C AB {P AB ) and C B v (Pbv) are defined by reference 
14 to similar CSF/vascular compliances as: 

C A b(Pab) = C AF (P AB ) Cab = 0.0209523 ml/mmHg 

C bv (Pbv) = C FV {P BV ) Cbv = 0.688845 ml/mmHg. 

[0055] When calculating extra-cranial compliances it is advantageous to first determine the 
total compliance of each compartment. This quantity will be denoted by C, with the subscript 
indicating the compartment. This type of compliance can be described by the relation 

Total Compartmental Compliance = Casein Compa rtmental Volume 

Increase in Compartmental Pressure 

and may be experimentally determined by measuring the pressure change that is induced by 
an estimated volume change and then taking the inverse of this ratio. Unfortunately, clinical 
data does not exist for determination of all of the compliances required for the current model. 
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However, there is information in the literature regarding volumes and distensibilities. These 
are related to compliance by 



Compliance = Distensibility ■ Volume. 



(79) 



[0056] Therefore, total compartmental compliance may be calculated from total 
compartmental volume (denoted Vj) and compartmental distensibility (denoted £>,-) through 
equation (79). However, before calculating compartmental volumes and distensibilities it is 
useful to describe how the resulting total compartmental compliances will be allocated to the 
local inter-compartmental compliances. 

[0057] Estimations of the local compliances are made by requiring that the sum of the 
relevant local compliances must equal the total compliance being apportioned. With this 
restriction, the model's predicted compartmental pressure changes induced by known volume 
changes should agree with experimental data. Based on this principle, the apportionment of 
total compliances among local compliances is described by 

Central, Local Compliances: 

C IY = Max[c } -C IT ,o] ml/mmHg (80) 
C YO = Max[c 0 -C TO ,o] ml/mmHg (81) 



C YM = Max[c Y - (C /Y +C ra ),0] ml/mmHg 



(82) 



Lower, Local Compliances: 



Czg = Cz ml/mmHg 



(83) 



Ccx = C,x ml/mmHg 



(84) 



C G m= Max[C G - (C ZG + C GX ), 0] ml/mmHg 



(85) 
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[0058] Here, the Max[ , ] function is used to ensure that no negative compliances occur in 
the model. For the particular choices of total compliances here, this function is not needed, 
but it is included above as it may be required if total tissue compliances are calculated to be 
much smaller. 

[0059] Total blood volume (V biood ) is estimated at 5000 ml and 5600 ml, where blood 
plasma makes up about 60% of this. The relative volumes of systemic veins, arteries and 
capillaries are: 



V sa = .16 Vbiood systemic artery volume (86) 

V$v 

.64 Vbiood systemic vein volume (87) 

V sc = .04 Vbiood systemic capillary volume (88) 

V pu i m = .09 V b i ood pulmonary system volume (89) 



V heart — -07 V b iood heart-blood volume (90) 

[0060] These values are based on estimates that the systemic arteries comprise 13% of the 
blood volume while the systemic arterioles and capillaries comprise 7%. In the present 
model, arteriole volume is lumped with the arteries, and hence 3% of the arteriole/capillary 
space is shifted here into the artery compartment resulting in the percentages above. It should 
also be noted that these percentages imply that systemic venous volume is four times that of 
systemic artery volume. 

[0061] As a next step, arterial and venous volumes must be allocated between the 
intracranial, lower, and central compartments. Estimates for this distribution are given 
below. It should be emphasized that these are percentages of the systemic blood and exclude 
the blood volumes in the heart and pulmonary system. 

pVcran = 07 percentage of systemic blood, intracranial region (91) 
pVi ower = .40 percentage of systemic blood, lower region (92) 
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pV central - -53 percentage of systemic blood, central region 



(93) 



[0062] The above percentages lead to the following volumes for the extra-cranial, vascular 
compartments: 

V z = pViower V sa = lower artery volume (94) 

Vi - pV r central V sa = central artery volume (95) 

V x =pVi oyver V sv = lower vein volume (96) 

V 0 = pV 'central Vsv = central vein volume (97) 

V D = pV !ower V sc = lower capillary volume (98) 



T r t r 



yj = py central Vsc = central capillary volume. (99) 

[0063] The volumes of the tissue and rest of body compartments are now calculated based 
on interstitial fluid volume and intracellular fluid volume. Total interstitial fluid volume is 
estimated in the prior art to be 12 L, while total intracellular fluid volume is estimated to be 
28 L. Considering that approximately 2 L of the intracellular fluid is contained in the blood 
cells, the following volume assignments can be made: 

V ^ = 12000 ml = interstitial fluid volume (1 00) 

V intra = 26000 ml = intracellular fluid volume (101) 

[0064] With these volume assignments and the percentages in equations (92) - (93), 
estimates for lower tissue and rest-of-body volumes become: 

V G = p Vio W er ( V ^er + Vintra) = lower tissue volume (102) 

Vy ■ — p V centra i ( V inter + V intra) + Vpulm rest of body volume (103) 
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[0065] Estimations for the pressure-strain modulus are given by Milnor and Nichols and 
O'Rourke for various branches of the systemic arteries. From this data, the distensibility of 
the various vessels can be calculated. Of interest here are the calculated distensibilities for 
the ascending aorta and the femoral artery, represented in the model by D f and D z . Milnor's 
data suggests that D f is 0.0036225 mmHg" 1 and D z is 0.00101 mmHg" 1 while Nichols and 
O'Rourke data suggests D } is 0.00320141 mmHg -1 and D z is 0.00237912 mmHg" 1 . 
Averaging these values results in the assignments 

Df = 0.00341 196 (104) 

D z = 0.00169456. (105) 

[0066] Since it is known that the central veins are eight times as distensible as central 
arteries, Do is calculated accordingly by 

Do=8D f . (106) 

[0067] There are some suggestions in the prior art that the upper venous distensibility may 
be larger that of the lower leg, but these results are inconclusive. Therefore, lower venous 
distensibility is set to 

D X =D 0 . (107) 

[0068] Consistent with prior art findings, the capillary compartments C, J, and D are 
considered non-deformable. 

[0069] The distensibility of the lower tissue compartment is assigned a weighted average of 
the interfacing vein and artery distensibilities. Since the veins have four times the volume of 
the arteries, these weights are used in 



Dq = ° z + 4 ° x = 0.0221754 (ml/mmHg)/ml (108) 
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[0070] Finally, the extensibility of the rest-of-body compartment (Y) is also a weighted 
average based on its composition: 



P V central inter + ^ intra + 



(A 



-V 



D v = 



pulm 



J 



D 0 + 



fl 



— V 



^5 



palm 



D 



J 



(109) 



= 0.0221838 (ml/mmHg)/ml 



[0071] The distensibilities of the extra-cranial compartments have now been determined. 
Combined with the volumes of the previous section, each extra-cranial compartment now has 
an associated total compliance equal to the product of its associated volume and 
distensibility. 

[0072] A lymphatic system is also included in the present model. This system allows 
interstitial fluid in the tissue compartments of the central and lower regions of the body to 
flow directly into the central venous compartment. Lymphatic flows are thus able to balance 
filtration in these regions and establish an equilibrium state. Lymphatic flow is not 
introduced in the intracranial region as the brain appears to lack lymphatic vessels. However, 
drainage of interstitial fluid via the Virchow-Robins Spaces is accommodated by pathways 
between the Brain and CSF compartments. 

[0073] The lymphatic system is represented in the present model through the flows Qyo and 
Qgo- These flows provide pathways through which interstitial fluid can move directly from 
the interstitium to the central venous system. Unlike the flows between adjacent 
compartments, which are driven by pressure differences, the lymphatic flow is governed 
almost exclusively by the interstitial fluid pressure. This flow increases 20 fold at interstitial 
fluid pressures near 4 mmHg from the corresponding flow at an interstitial pressure of -6 
mmHg. Below this interstitial fluid pressure value, lymphatic flow ceases. The relative 
lymphatic curve is modeled by the logistic expression 



Plymph(x) = M r(x+6y where M= 20 and r = 1 (110) 

F \ + (M -\)e~ r(x+6) 
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[0074] This function, depicted in FIG. 4, defines the relative lymph flow in terms of 
interstitial pressure. The large dot in FIG. 4 denotes the mean pressure of interstitial fluid in 
the model. Conversion of this relationship to the actual lymphatic flows Qgo and Q Y o is 
accomplished by defining 



[0075] In this manner, the mean flow is maintained at the mean pressure and can increase 
or decrease accordingly. Based on equation (110), Qgo can increase from a mean flow of 
about 0.24 ml/min to 4.7 ml/min while Q Y o can increase from 1.76 ml/min to 35 ml/min. 

[0076] Blood supply and CSF production is well regulated in the human intracranial 
system. This auto-regulation is achieved by the dilation and constriction of the blood vessels 
induced by factors such as carbon dioxide concentration, hydrogen ion concentration, oxygen 
concentration and intracranial pressures. The model incorporates such auto-regulation in the 
form of pressure sensitive fluidities depicting either vessel constriction in the form of a 
reduced fluidity or dilation in form of increased fluidity. 

[0077] In all, the present model contains 16 compliance parameters, of which eight are 
variable and depend on pressure differences. There are also 23 resistance parameters. Four 
intracranial resistances are pressure-dependent and autoregulate flow to both the cerebral 
capillary bed and the choroid plexus as well as the production of CSF in the choroid plexus. 
CSF production is known to be nearly constant over a wide range of pressures. However, the 
production 35 of CSF decreases when perfusion pressure decreases to low values, and a 
mechanism that ramps down CSF production when perfusion pressure falls below 55 mmHg 
is also included in the model. 

[0078] Blood flow into the cerebral capillaries is auto-regulated by a pressure dependent 
fluidity Z AC defined by 



Qgo(0 



p iymph ( p g) Qgo 

P lymph ( P G) 



(111) 



QyoO) 



Pjymph ( P Y ) QYO 
P tymph (Py ) 



(112) 
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(113) 



where F(P A ) is a logistic expression defined by 




- with M= 1 .0000 1 and r = 0.25 . 



(114) 



[0079] This implies that Q AC = Zac (Pa-Pc) = F (Pa) and cerebral blood flow is thus 
determined almost entirely by intracranial artery pressure and does not diminish significantly 
until this pressure drops below 50 mmHg. The relationship between Q AC and P A is displayed 
in FIG. 5. 

[0080] Blood flow into the choroid plexus is regulated by a pressure-difference dependent 
fluidity Z AP defined by 



where the multiplier G(P perf ) in equation (115) is a function of the perfusion pressure P P er/= 
P A -P B that is unity when P perf exceeds 55 mmHg and then falls linearly to zero below this 
value. This implies that 



[0081] The multiplier G has been included in equation (115) to model the fact that CSF 
production in the choroid plexus scales with blood flow, remains nearly constant over a wide 
range of pressures, and is decreased as a function of the magnitude of the perfusion pressure 
when P pe rf falls 20 below 50 to 60 mmHg. From equation (116), blood flow into the choroid 
plexus remains constant unless perfusion pressure falls to a low level. 




(115) 



Qap=Zap ' (Pa-Pp)^Qap • G(p perf ). 



(116) 
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[0082] Over the wide range of normal pressures, the production of CSF is auto-regulated at 
the venous side of the choroid plexus capillaries by the pressure dependent fluidity Z PV 
defined by 

„ QAP ' G { P perf)- Z PF ' { P P ~ P f) ni « 
^PV ~ ■ U 1 ') 



[0083] When P per f> 55, this expression for Z P v will maintain a constant pressure difference 
between the choroid plexus and the ventricular CSF. Substituting Z PV with G = 1 into the 
governing equation for the choroid plexus, equation (15) reduces to 



CpF ^£^ = 0. (118) 
dt 



[0084] Since the compliance C PF must be non-zero to account for the known ability of the 
Choroid Plexus to transmit pressure pulsations to the ventricular CSF, the governing equation 
for the choroid plexus compartment becomes simply 



—^ = 0. (119) 
dt 



[0085] This implies a constant pressure difference between the choroid plexus and 
ventricular CSF is maintained by equation (1 17) for P per f> 55 mmHg. Therefore 

Q PF =Z PF • (P P - P F )= Z PF (P P -P F )=Q PF . (120) 



[0086] Since for pressures in the normal range, CSF production in the choroid plexus is 
proportional to P PF , constant CSF production from the choroid plexus is thus achieved. 
Equation (119) also eliminates the need to estimate C P p in this model as occurrences of this 
parameter are always multiplied by either dP PF /dt or dPpp/dt. 
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[0087] Using the above results, it can now be demonstrated how Z PV autoregulates CSF 
production for P per f> 55 mmHg. Substituting Q PF for Z PF • (P P - P F ) in equation (1 1 7) and 
noting that Q AP -Q PF = Q PV results in the equalities: 



[0088] The last term in this expression reveals the relationship between Z PV and P F . 
Physiologically, Z PV should decrease with increasing CSF pressure (P F ) causing an increase 
in choroid plexus pressure and maintaining constant CSF production across Z PF . It is quite 
clear from the last equality in equation (121) that this is indeed the case. These 
autoregulatory mechanisms may be substituted into (13) to (16) and (18) to obtain the 
governing equations for compartments A, C, P, V, and F. 

[0089] The present model also includes a group of regulatory mechanisms associated with 
the sympathetic nervous system (SNS). Two variable resistances in the central and lower 
regions provide for SNS regulation of arterial pressure through constriction of the arterioles. 
Included in these variable resistances is the dramatic "last ditch stand" triggered by the SNS 
when arterial pressure falls to a level where cerebral blood flow is affected. The far less 
extreme SNS regulation of arterial pressure through a constriction of the large vascular 
vessels in the central body is also represented in the model. Active interfaces are placed 
between the central rest-of-body compartment and the central artery and vein compartments. 
When arterial pressure falls, forcing terms in the governing equations for compartments Y, I, 
and O force the volume cups at the active Y-I and Y-O interfaces into the vascular 
compartments, providing the regulatory constriction. An additional SNS mechanism in the 
model regulates central arterial pressure by increasing the number of heartbeats per minute if 
arterial pressure falls. 

[0090] The sympathetic nervous system (SNS) is associated with reflex mechanisms that 
act to maintain levels of arterial pressure and cardiac output if arterial pressure falls. This 
section deals with the portions of the model that represent the capacity of the SNS to rapidly 
increase arterial pressure by constricting both the arterioles and the large vascular vessels. 
The model's regulation of cardiac output by the SNS will be developed below. 



Z P y = 



Qpv 




(121) 



Pp-Py 



(P P -P F )+(P F -Py) Q PF /Z PF +{P F -P v y 
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[0091] A first SNS pressure mechanism involves constriction of the large vessels of the 
circulation in the central body. If arterial pressure falls, the large vessels of the circulation in 
the central body, especially the veins, strongly constrict to cause a rapid increase in arterial 
pressure. This mechanism is included in the current model by placing active interfaces 
between the central compartment Y and the central vascular compartments I and O. Forcing 
terms in the conservation equation for compartment Y now force the volume cups at the 
interfaces to push into the I and O compartments when arterial pressure diminishes, modeling 
the SNS large vessel constriction mechanism. 

[0092] The conservation equation in compartment Y is described by 



at dt dt dt 



where the second equality determines the equation. Here, the volume change dV YM fdt 
between compartment Y and the external environment M, as in equation (3), is simply 
proportional to the change in pressure difference dPy^/dt. However, dVyc/dt and dVy/dt 
involve both changes in the pressure differences dPyo/dt and dPy/dt, respectively, and a 
forcing term describing the active compliance between Y and O and between Y and I, 
respectively. In particular, 



and 



dV 



YO 



dt 



= C 



dP 



YO 



YO 



dt 



+ F Yo(* ) with F vo(t) = -24P'/ 



(123) 



dV 



YI 



dt 



= c 



dp 



YI 



YI 



dt 



+ F YI (t ) with F YI (t) = -8P7 



(124) 



[0093] The governing differential equation in compartment Y is now obtained by replacing 
dVyo/dt and dV YI /dt in equation (7) by the expressions (123) and (124). Similarly, for the 
governing equations in compartments O and I, dVyo/dt is replaced in equation (8) by 
expression (123) and dVy/dt is replaced in equation (5) by expression (124), respectively. 
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[0094] A second SNS pressure regulation mechanism involves constriction of the arterioles 
in most parts of the body (excluding the brain and heart) when central artery pressure drops. 
This causes an increase in arterial pressure through an increase in the total peripheral 
resistance. When the arterial pressure drop is severe enough to compromise blood flow to the 
brain, this regulatory response is extreme. 

[0095] To model this portion of the SNS regulatory response, two multipliers are defined 
for the artery-capillary fluidities Z/j and Z Z d- These variable resistances are of the form 



Z " = • SNSz x (Pj) • SNSz(Q AC ). (125) 



and 



Z ZD = » QZ % • SNSztfPj) ■ SNSz(Q AC ). (126) 



[0096] The first multiplier, 



SNSz } (Pj) = %-, (127) 

n 



is a function of central body artery pressure that increases resistance if Pj falls below its mean 
level. The second multiplier, which is a function of the cerebral blood flow Q AC , is defined 
by 



SNSz(Q AC ) = v where M= 1.1 and r = 0.02. (128) 

\ + (M-\)e r y Q ^-QAC) 



[0097] This multiplier models the last-ditch stand when cerebral blood flow is significantly 
reduced. SNSz remains close to unity (not activated) until arterial pressure drops to a level 
where Q AC is affected. At this point, SNSz drops sharply, dramatically increasing the 
resistances R/j and Rzd- 
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[0098] The major source of forcing in the present whole-body model comes from the heart 
and involves the cardiac output Qhi and cardiac uptake Qoh- Because it is known that all 
extra blood that flows into the heart is automatically pumped without delay into the aorta and 
flows again through the circulation, the cardiac output will be set equal to cardiac uptake, so 
that Qhi = Qoh- It is further known that the normal heart, functioning without any special 
stimulation, can pump an amount of venous return up to 2.5 times the normal venous return 
before the heart becomes the limiting factor. Therefore, a venous return function R is defined 
that incorporates all of the flow into the central venous compartment 

R = Qyo +Qgo + Zj 0 (Pj -P 0 ) + Z TO (Pt ~Po) + Z so (Ps -Po) + Z xo (P x ~Po\ (129) 

[0099] The cardiac uptake Qoh is now defined in terms of this venous return and two 
regulatory multipliers SNSo and OVP by 

Qoh = SNSo (P f j) • OVP(P Q ) ■ R (130) 

where P 7 denotes the time derivative of the arterial pressure function Pj . Since Qhi = Qoh, 
cardiac output is now based on the venous return through equation (130). 

[0100] The multiplier SNSo in equation (130) is associated with cardiac regulation by the 
sympathetic nervous system (SNS). It models an SNS mechanism that increases or decreases 
the heart rate in response to a change in arterial pressure. A decrease in arterial pressure can 
markedly increase heart activity, with the heart rate increasing to as high as 160 - 1 80 bpm 
from the normal value of 72 bpm. A linear model for the SNSo multiplier as a function of the 
instantaneous time derivative P'j is 

SNSo(P'i) = \ - i (131) 

10 

[0101] This relationship produces an increase in heart rate to about twice the normal level 
as the result of a pressure drop of 20 mmHg. Note that SNSo (0) = 1, so that this regulatory 
mechanism is only active when Pj is changing. 
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[0102] The OVP function in equation (130) insures that if venous pressure drops, then so 
does cardiac uptake. It is defined by 



OVP(P 0 ) = 



M 



where M= 2.5 and r — 0.5 . 



(132) 



\+{M-\y{ p o- p o) 




Since OVP\Pq) = 1, this regulatory mechanism is not active when central venous pressure 
remains at its mean value. 

[0103] Using the above cardiac forcing terms leads to a system of equations that describe 
mean pressures in the sense that the oscillations that occur about a mean pressure value are 
removed and the mean pressure may be considered a temporal average over one cardiac 
cycle. However, when trying to resolve circulatory pressure pulsations caused by cardiac 
output, a more instantaneous description of cardiac output is necessary. Again in this case, 
cardiac output is set equal to cardiac uptake. 

[0104] Referring now to FIG. 6, it is contemplated that mathematical model 10 of the 
present invention will operate in a computing environment including a computer system 20. 
Computer system 20 includes a computer central processing unit (CPU) 22, a computer 
memory 24, and input/output devices 26. Mathematical model 10 is typically implemented 
with or embodied in a computer program 28 which, when executed by computing resources 
within computer system 20, provide the functionality of the present invention. Typically, 
computer program 28 resides in computer memory 24 of an individual client computer 
system. Of course, computer program 28 may reside in the memory of a local or wide area 
network server or in the memory of an equipment computer processor. In one embodiment, 
computer programs 28 include a body compartment module 29, a flow equation module 30, 
and an equation solver module 31. Body compartment module 29 includes instructions for 
dividing the body into a plurality of compartments and a heart pump, each compartment 
representing a portion of the body. Flow equation module 30 includes instructions for 
deriving a plurality of differential flow equations, each of which correspond to one of the 
compartments. Equation solver module 31 includes instructions for solving the plurality of 
differential flow equations. The instructions within body compartment module 29, flow 
equation module 30, and equation solver module 31 are executed within computer programs 
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28 to simulate the pressure dynamics of the intracranial system. Input/output devices 26 
typically include a storage device 36, such as a hard disk drive, a keyboard 38, a pointing 
device 40, i.e., a mouse, and a display device 42, such as a monitor. It is contemplated that 
data may be manually input to mathematical model 10 via input/output devices 26 or 
automatically input via patient monitoring equipment or the like. 

[0105] To validate the present whole-body model for intracranial pressure dynamics, 
several types of simulations assuming normal physiology were carried out and the results 
compared to clinical observations. In the first type of simulation, a pulsatile form for cardiac 
output was introduced as forcing so that the model's predicted pressure responses could be 
compared to clinically measured pressure pulsations in healthy humans. The pulsatile cardiac 
output function with parameter values n = 1 1, cp = .27882, and a = 5.958 was used as the 
cardiac forcing function (Qui) in the model's governing equations. This function is depicted 
in FIG. 7 A and represents the cardiac output over two cycles. A mean value for central artery 
pressure of 96 mmHg was prescribed as an initial condition for the first simulation. In the 
second type of simulation, a constant flow infusion test was used to validate the model's 
representations for variable intracranial compliances and intracranial fluid dynamics. In these 
simulations, the governing equation for the extra-ventricular CSF compartment (T) was 
augmented by adding a constant infusion term to model the clinical infusion of mock CSF 
into the lower lumbar space. 

[0106] With the lower body region explicitly represented by separate compartments in the 
present model, the effect of orthostatic forces on cerebral blood flow as the body changes 
position can be considered. Two additional types of simulations were now run to validate the 
modeling of the regulatory mechanisms associated with the systemic nervous system. The 
normal physiology value for the resistance R X o was increased twofold and R IZ was decreased 
by one half to simulate a change in body position from lying down to standing up, and the 
behavior of the cerebral blood flow Q AC was then examined. These resistance changes were 
made instantaneously and the modeled effect on cerebral blood flow was determined for the 
next 30 seconds. In the first simulation of this type, the model equations included all of the 
SNS reflexes described above. In the second simulation, the SNS terms were removed from 
the governing equations. 
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[0107] To examine the potential predictive capabilities of the model in pathological 
conditions, simulations were run where all cardiac forcing in the model was suddenly stopped 
to simulate cardiac arrest, i.e., Qm and Qoh were suddenly set equal to zero. A blood volume 
of 5600 ml was assumed in these simulations, and the response of pressures in the model's 
circulatory compartments was determined. This behavior, and the predicted final circulatory 
compartmental pressure values, were then compared to clinical results associated with cardiac 
arrest and the filling pressure of the circulation. 

[0108] As a second example of pathology, simulations of hemorrhagic shock were carried 
out. Hemorrhage was modeled by the inclusion of an outflow path, denoted Q X m, from lower 
venous compartment X into ambient atmosphere compartment M. This flow was calculated 
so as to achieve a 45% loss in blood volume at the end of the simulation. The percent 
changes in central artery pressure, cardiac output, and cerebral blood flow were then 
calculated with respect to percent blood loss. 

[0109] In all of these simulations, the model's system of differential equations was solved 
numerically using a typical symbolic mathematical software package employing maximum 
accuracy settings. 

[0110] The usual systolic and diastolic values of pulsatile central artery pressure are 
perhaps the best-known values in human physiology. FIG. 7B shows the behavior of central 
artery pressure predicted by the present model in response to the pulsatile cardiac forcing 
developed in the prior art. The model's predicted response to a mean value of 96 mmHg, 
prescribed as an initial condition, reproduces a peak systolic pressure of 120 mmHg and a 
diastolic pressure of 80 mmHg. FIG. 7B also depicts the intracranial artery pressure 
response. Here, the predicted systolic and diastolic pressures are "100 over 65." These values 
are typical of those in the small arteries. 

[0111] Because even the most complex mathematical model must be based on assumptions 
and simplifications of actual physiology, model validation is an essential step in the 
development process. In the present case, after calibration of parameters associated with 
healthy human physiology, the model was used in two types of simulations and the results 
compared to physical data. In the first type of simulation, the response of compartmental 
pressures to the realistic pulsatile cardiac output given in FIG. 7 A was determined. As shown 
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in FIG, 7B the predicted response of central arterial pressure is the typical "120 over 80" 
blood pressure readings expected in healthy humans. It is worth noting that a mild incisura 
(dicrotic notch) is discernable in this pressure curve. This notch is typical in central artery 
pressure readings, although is it usually more pronounced than in FIG. 7B. In particular, the 
present model cannot capture reflected waves which have a tendency to cause a second 
(usually higher) systolic peak and hence a more prominent inflection point. 

[0112] The response of the intracranial arteries is also shown in FIG. 7B. This response 
agrees well with prior art estimates for pressure pulses in the small arteries. The pressure 
responses of other intracranial compartments were also within expected ranges. 

[0113] Constant flow infusion tests were also simulated using the present model. In these 
experiments, mock CSF was infused at a constant rate into the lower lumbar space. The 
pressure of this space was then measured and associated with a calculated total CSF volume 
change. This resulted in determination of a curve known as the global pressure-volume 
relation. The slope of this curve describes the elastance of the entire CSF space, including 
extracranial portions. The inverse of the elastance is the more well-know compliance. 

[0114] The typical clinical global CSF pressure-volume relation, except at extreme 
pressures, is an S-shaped curve of logistic type. It has a lower pressure plateau near resting 
pressure, defined as the pressure where CSF production in the system is just balanced by CSF 
absorption by the venous system. This region of small slope (large compliance) is due to the 
ability of the system to easily accommodate increases in the volume of the CSF space at these 
relatively low pressures through the compression of the venous system (excluding the venus- 
sinus veins). As additional CSF volume is added to the system and pressures increase, this 
capacity for adjustment diminishes as there will be less venous blood available to eject to 
further compress the veins. Thus, with increasing pressures, the pressure-volume curve 
steepens indicating a reduction in the compliance of the system. The slope of the pressure- 
volume relationship continues to increase for larger infusion volumes until the resulting CSF 
pressures are high enough that the intracranial arteries can begin to be compressed. At this 
point, some additional compliance enters the system. There is a point of inflection in the 
curve followed by a region of decreasing slope that leads to an upper pressure plateau at the 
diastolic pressure of the intracranial arteries. Once additional volume increases beyond this 
point increase CSF pressures to the systolic pressure of the intracranial arteries, there are no 
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additional mechanisms available to buffer additional volume increases, and the compliance of 
the CSF system falls to zero. 

[0115] The pressure-volume relation obtained by using the present model to simulate CSF 
infusion tests is given in FIG, 8A. The pressure in this figure is the pressure P T of the 
extraventricular CSF compartment since this compartment contains the lower lumber portion 
of CSF space where clinical pressure recordings were made. The volume change in FIG. 8A 
represents the net change of both ventricular and extra- ventricular CSF volumes combined, as 
this is the quantity calculated in the experiments. The logistic-like shape of the predicted 
pressure-volume relationship is maintained until CSF pressures reach high values comparable 
to the systolic blood pressure. At this point, the curve turns upward and compliance rapidly 
decreases. 



[0116] The same infusion simulations that led to the pressure-volume relationship in FIG. 
8A may also be used to calculate a relationship between CSF pressure increases and CSF 
absorption. The slope of this relationship is known as the conductance of CSF outflow and is 
denoted by C out . In model simulations, this value is easily calculated by 



r m (A Z TS {P T {t)-P s (/ )) + Z TO (P T (t)-P Q (Q) - [Qrs + Qtq ) 

c ou,vn- p T {t)-p T K } 



where the superscript m denotes the model calculation and t denotes the time in minutes from 
the start of the simulated infusion. FIG. 5B illustrates the predicted conductance of CSF 
outflow throughout ten minutes of an infusion simulation. 

[0117] When C out has been has been studied experimentally, a linear relationship between 
CSF pressure increase and CSF absorption is observed. In prior art experiments, a mean 
value for C out of 0.1 1 (ml/min)/mmHg is given for a sample of eight healthy volunteers, and it 
is stated that a value greater than 0.10 is probably normal. The values of the conductance of 
CSF outflow calculated from the present simulations are shown in FIG. 8B. These values 
change with time, but stay within .004 units of the mean value of 0.1 1 (ml/min)/mmHg 
observed in the prior art. They are also greater than 0.10 as previously suggested. 

Furthermore, the calculated temporal variation of C^ ut is sufficiently small that that the 
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relationship between CSF pressure increase and CSF absorption might easily be categorized 
as linear on the basis of clinical data. 

[0118] In the current model calibrations, the pressures, flows, and hence resistances are 
determined from data associated with the body being in the horizontal (supine) position. 
Gravitational influences on the circulation will change when the body changes position. 
Upon standing up, the blood flow into the lower body is aided by gravity, but the venous 
return from the lower body is hindered by it. As the lower body in the present model is 

■ 

represented by separate compartments, the effect of a positional change can be considered. 
To simulate the gravity-induced changes associated with standing up, the resistance into the 
lower arteries (R IZ ) may be decreased by one half while the resistance from the lower veins 
(Rxo) is doubled. This results in pressure and volume increases of the lower arteries and 
veins. In the current simulations, these resistance changes were made instantaneously and the 
modeled effect on cerebral blood flow was determined for the next 30 seconds. In the first 
simulation, all sympathetic nervous system (SNS) reflexes are intact and in the second 39 
these are removed. As can be seen in FIG. 9, with all SNS reflexes activated, cerebral blood 
flow is maintained at 100% of its original value. 

[0119] Referring now to FIG. 9, the results of modifying the resistances between the central 
and lower body in the model to simulate a change in body position from lying down to 
standing up are illustrated. The solid curve in FIG. 9 indicates the response of cerebral blood 
flow to this change in position with the SNS reflexes included in the model equations. The 
simulation predicts that with all SNS reflexes activated, cerebral blood flow is maintained at 
100 % by the model through the change in body position. The dashed curve in FIG. 9 shows 
the predicted behavior of cerebral blood flow in the absence of regulation by the SNS. 

[0120] When the SNS reflexes are removed from the model equations, cerebral blood flow 
drops to about 94% of its original value due to the positional change. This predicted decrease 
indicates the important role that the SNS regulatory mechanisms play in the maintenance of 
cerebral blood flow in the model. The decrease also confirms that loss of the SNS reflexes 
can be a factor in orthostatic intolerance. Syncope, or fainting, is a common event with many 
different causes. Syncope can be due to sudden vasodilatation (vasodepressor or "vasovagal" 
syncope) as well as postural hypotension, in which the normal vasoconstrictive reflex 
response to a transiently decreased cardiac output on standing is not sufficiently active. 
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Postural hypotension can occur as the result of drugs, venous disease, sympathectomy, 
hypovolemia, peripheral neuropathy, in addition to degeneration of the sympathetic nervous 
system (primary autonomic insufficiency, or idiopathic orthostatic hypotension). Common to 
all causes of syncope is a decrease in cerebral blood flow to under 30 ml per 100 g brain 
tissue per minute from the usual 50-55 ml, about 55% of normal. This is well above the 
threshold for loss of electrical function, at 30% of normal, but is enough to cause a transient 
loss of consciousness. The calculated percentage decrease- in cerebral blood flow in the 
current simulation with the SNS terms in the model equations omitted does not approach the 
levels associated with fainting due to a rapid change in position as the result of standing up 
too quickly, even with inactive SNS reflexes. However, cerebrovascular autoregulation in 
the current simulation remains uncompromised and will act to maintain cerebral blood flow 
despite sudden decreases in arterial pressure. 

[0121] To demonstrate the potential predictive capabilities of the model in pathology, two 
situations were simulated. In the first of these, a simplified representation of cardiac arrest 
was created by suddenly terminating all cardiac forcing, i.e., setting Qhi = Qoh = 0, at a 
specified time in the course of the simulation. In prior research by Gutyon, it is noted that 
"When heart pumping is stopped ... the flow of blood everywhere in the circulation ceases a 
few seconds later. Without blood flow, the pressures everywhere in the circulation become 
equal after a minute or so. This equilibrated pressure level is called the mean circulatory 
filling pressure ... at a volume 40 of 5000 ml, the filling pressure is the normal value of 7 
mmHg." Predicted results from the simulation were consistent with this statement. When all 
cardiac forcing was stopped in the simulation, all circulatory pressures tend to 7.5 mmHg 
after about 45 seconds. The small difference in the predicted and cited filling pressures may 
be partially due to the fact that total blood volume in the simulation was taken to be 5600 ml 
rather than the 5000 ml. 

[0122] Referring now to FIG. 10, the second pathological situation considered as a test of 
the model's capabilities involves the effects of shock caused by hypovolemia. In these 
simulations, a flow term Q X m was introduced into the equation for compartment X to model a 
hemorrhage from the lower body. Prior art research by Guyton gives a clinically-derived 
graph that, as in FIG. 10, depicts the percentage drop in central artery pressure (Pj) and 
cardiac output (Q H i) with respect to percent blood loss. Consistent with these clinical results, 
the modeled response shows very stable pressure and cardiac output over the first 10% blood 
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loss, and the relative arterial pressure stays above relative output during the entire course of 
the hemorrhage. At 30% blood loss, the relative drops in pressure and cardiac output are also 
in close agreement with the clinical results. During the course of the simulated hemorrhage, 
the central venous compartment contracts to less than 50% of its original volume, indicating 
that the active interface (SNSc) in the model is causing the central veins to strongly constrict 
in order to maintain blood flow back to the heart. Beyond 30% blood loss, a second arterial 
plateau is noticed in FIG. 10. This important feature, which is also noted in the clinical 
results, is due to the "last-ditch stand" reflex to maintain blood supply to the brain provided 
by the regulatory multiplier SNSz (Q A c) in equation (130). 

[0123] Referring now to FIG. 1 1, the predicted percentage drop in cerebral blood flow 
(Qac) with respect to percent blood loss during the simulated hemorrhage is illustrated. The 
solid curve represents results when both SNSz mechanisms are intact, and the dashed curve 
represents results when both SNSz mechanisms are disabled. Both simulations predict that a 
25% blood loss can be tolerated without a significant reduction in cerebral blood flow. 
Above this degree of blood loss, with the SNS mechanisms intact, the blood supply to the 
brain remains above 75% of its original value. However, without these SNS regulatory 
mechanisms included in the model's equations, cerebral blood flow drops quickly to 20% its 
original value. 

[0124] The role of the SNS regulatory mechanisms in maintaining adequate central artery 
pressure and cardiac output has been noted above. However, the need to include the arteriole 
constrictive reflex SNSz in a model of intracranial pressure dynamics when simulating 
pathology is also dramatically illustrated by the behavior of cerebral blood flow as a function 
of blood loss in the hemorrhage simulations. Due to the intracranial autoregulatory 
mechanism Z AC in the model, a 25% blood loss can be tolerated without significant reduction 
in cerebral blood flow. Beyond this point, arteriole constriction is essential to further 
maintain adequate levels of cerebral blood flow. As indicated in FIG. 1 1, predicted cerebral 
blood flow Q AC remains above 80% of its original value from a 25% blood loss through a 
40% blood loss. By contrast, in an analogous simulation where the arteriole constrictive 
response was disabled, Q AC falls to approximately 55% of its original value when a 40% 
blood loss has occurred. Between a 40% and 45% blood loss, Q AC remains above 75% of its 
original value with the arteriole constrictive response intact, but falls rapidly from 55% to 
approximately 20% of its original value with the arteriole constrictive response disabled. 
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[0125] Most attempts to study intracranial pressure using lumped-parameter models have 
adopted the classical "Kellie-Monro Doctrine," which considers the intracranial space to be a 
closed system that is confined within the nearly-rigid skull, conserves mass, and has equal 
inflow and outflow. The present work revokes this Doctrine and develops a mathematical 
model for the dynamics of intracranial pressures, volumes, and flows that embeds the 
intracranial system in extensive whole-body physiology. The new model consistently 
introduces compartments representing the tissues and vasculature of the extradural portions 
of the body, including both the thoracic region and the lower extremities. In addition to 
vascular connections, a spinal subarachnoid cerebrospinal fluid (CSF) compartment bridges 
intracranial and extracranial physiology allowing explicit buffering of intracranial pressure 
fluctuations by the spinal theca. The model contains cerebrovascular autoregulation, 
regulation of systemic vascular pressures by the sympathetic nervous system, regulation of 
CSF production in the choroid plexus, a lymphatic system, colloid osmotic pressure effects, 
and realistic descriptions of cardiac output. 

[0126] Validation tests show that the model's predictions are in agreement with 
experimental data for normal physiology. Additional simulations indicate that the present 
whole-body model appears to have significant potential predictive capabilities in situations 
involving pathological conditions. 

[0127] One embodiment of the present invention has been disclosed and discussed herein, 
although it should be understood that the present invention is not limited to this (or any other) 
particular embodiment. On the contrary, the present invention is intended to cover all 
alternatives, modifications and equivalents that may be included within the spirit and scope of 
the appended claims. 
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